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Abstract 



Results for the AA^ and UN interactions obtained at next-to-leading order in chiral effective field theory are reported. 
^^ . At the order considered there are contributions from one- and two-pseudoscalar-meson exchange diagrams and from 

^^ ' four-baryon contact terms without and with two derivatives. SU(3) flavor symmetry is imposed for constructing the 

0> . hyperon-nucleon interaction while the expUcit SU(3) symmetry breaking by the physical masses of the pseudoscalar 

mesons (tt, K, rj) is taken into account. An excellent description of the hyperon-nucleon system can be achieved at 

next-to-leading order It is on the same level of quality as the one obtained by the most advanced phenomenological 

hyperon-nucleon interaction models. 
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1. Introduction 

'^ ' While there is a steady interest in physics involving baryons with strangeness and a corresponding increase of 

-Y-\ ■ empirical information ever since the discovery of the A-hyperon many decades ago, the present times seem to be 

f^ ' particularly rewarding. First, at new experimental facilities like J-PARC in Japan or FAIR in Germany a significant 

lO ■ amount of beam time will be devoted to strangeness physics research. The proposed experiments encompass accurate 
measurements of level spectra and decay properties of strangeness S - -I and S - -2 hypernuclei fl, 2j but also of 

f^ . elementary cross sections for Y.^p scattering [3]. Information on Kp scattering, specifically on the scattering lengths, 

CO ' might emerge from ongoing studies of the final-state interaction in production reactions like pp — > K^Ap |4] and 

Parallel to this development, techniques for dealing with few- and many-body systems have reached a high degree 
of sophistication ia, tZl Isl |9l [lOl [ijj, llS, [131]. Some of these allow one to consider nuclei with much more than four 



Lj . oi sopinsucaiiou iimi/L ioli^ ihaii iliizi.ii jii . ooiiie oi uiese aiiow one lo coiisiuci iiucici wiui iiiucii iiiuic man luui 

i*^ ' nucleons, the limit for standard few-body calculations with the Faddeev- Yakubovsky theory [7] . Of particular interest 

C^ . in this context are nuclear lattice simulations as they offer a new many -body technique directly tailored to the effective 

field theory description of baryon-baryon interactions, as high-lighted recently by the first ever ab initio calculation of 
the Hoyle state in the spectrum of '^C [14]. Thus, it seems to be feasible to perform similar calculations of hypernuclei 
too, with comparable accuracy as those for ordinary nuclei, which would open a completely new testing ground for the 
hyperon-nucleon ( YN) interaction. Though few -body calculations of hypernuclei can be akeady found in the literature 



015111611 . for the latter aspect it would be desirable to employ techniques that allow one to use directly the elementary 
YN interaction (i.e. without any approximation) and, in particular, to include the important A-2 conversion. Only 
then one can connect the properties of the hypernuclei unambiguously with those of the underlying AN (and UN) 
interaction. 

Finally, and on a different frontier, lattice QCD calculations have matured to a certain degree, as documented in 



recent review articles Ill7[ 11811 . and are coming closer to a level where they can provide additional constraints on the 



baryon-baryon interactions in the strangeness sector |119 1 



Preprint submitted to Elsevier April 22, 2013 



To keep up with these developments we present here a study of the YN interaction performed at next-to-leading 
order (NLO) in chiral effective field theory (EFT). It builds upon and extends a previous investigation by the Bonn- 
Jiilich group carried out at leading order (LO) [.20J . Using chiral EFT for the YN interaction is prompted by the great 
success that this scheme has met in the application to the nucleon-nucleon {NN) interaction. Indeed, proposed by 



Weinberg 1121112211 more than two decades ago, chiral EFT has turned out to be a rather powerful tool for the derivation 
of nuclear forces. Its most salient feature is that there is an underlying power counting which allows one to improve 
calculations systematically by going to higher orders in a perturbative expansion. In addition, it is possible to derive 
two- and three-nucleon forces as well as external current operators in a consistent way. We are now at a stage that the 
latter aspect will also be important for realistic studies ofhypernuclear interactions. In the past there have been already 



discussions on the role of a three-body YN interaction 0231 12411 in hypemuclei and specifically for the properties of 



neutron star matter Il25[ 12611 



As the most recent applications demonstrate, the nucleon-nucleon interaction can be described accurately within 



the chiral EFT approach 11271 12811 . In line with the original suggestion of Weinberg, the power counting is applied to 
the A^A^ potential rather than to the reaction amplitude. The latter is obtained from solving a regularized Lippmann- 
Schwinger equation for the derived interaction potential. The chiral A^A^ potential contains pion-exchanges and a 
series of contact interactions with an increasing number of derivatives. The latter represent the short-range part of the 
A^A^ force and are parametrized by low-energy constants (LECs), that need to be fixed by a fit to data. For reviews we 



refer the reader to the recent Refs. [l29[|30 l 



In our study of the FA^ interaction, we follow the scheme that has been applied by Epelbaum et al. 0281 l3ll 13211 
to the A^A^ interaction. For investigations of the FA^ interaction based on other schemes see 113311341 1. Still, there are 
some essential differences between the AA^, EA^ systems and the A^A^ case that have an influence on how one proceeds 
in the application of chiral EFT in practice. First and foremost, there is no phase-shift analysis for the S - -I sector 
and, therefore, we have to fix the LECs by a direct fit to data rather than by a fit to individual partial waves as it is 
done in the A^A^ case. Secondly, the amount of YN data is rather limited. Indeed, there are basically only integrated 
cross sections, often with large uncertainties. Thus, we follow here the practice of previous investigations of the YN 



interaction, notably those performed in the meson-exchange picture 11351 1361 1371 13811 . and impose constraints from 
SU(3) flavor symmetry in order to reduce the number of free parameters. In particular, all the baryon-baryon-meson 
coupling constants are fixed from SU(3) symmetry and the symmetry is also exploited to derive relations between 
the various LECs. In the actual calculation the SU(3) symmetry is broken, however, by the mass differences between 
the Goldstone bosons {n, K, tj) and between the baryons. For these masses we use the known physical values. In 
any case, we want to stress that we consider the imposed SU(3) symmetry primarily as a working hypothesis and not 
as a rigorous constraint. Future data with higher precision will possibly demand to depart from SU(3) symmetry in 
some way. In that sense our present investigation certainly has primarily an exploratory character At the moment 
we are able to describe the available AN and ILN data consistently without any explicit SU(3) breaking in the contact 
interactions as will be demonstrated below. A simultaneous description of the A^A^ interaction with contact terms that 
fulfil SU(3) symmetry turned out, however, to be not possible. 

As mentioned above, in order to obtain the reaction amplitude from the interaction potential derived within chiral 
EFT, one has to solve a regularized Lippmann-Schwinger equation. The question how this regularization should be 
performed is an open issue and is still controversially discussed in the literature, see, e.g. [39, 40, A\j 42i- In the 
present work, we refrain from touching this certainly very important question. Rather we follow closely the procedure 
adopted by Epelbaum et al. 128113211 and others [27], in their study of the A^A^ interaction and introduce a momentum- 
dependent exponential regulator function into the scattering equation. 

The present paper is structured as follows: In Sect. 2, a review of the chiral EFT approach is given with special 
emphasis on the imposed SU(3) symmetry. In particular, the structure of the contact interactions at LO and NLO is 
specified and the expression for the one-meson exchange contributions are reproduced. A detailed description of the 
two-boson exchange potential that arises at NLO is presented in the appendix. The strategy followed in the fit to the 
data is outlined in Sect. 3. In Sect. 4 results for the AN and I,N interactions obtained at NLO are discussed and com- 



pared to available experimental information. Results of our LO calculation and of the Jiilich '04 FA^ interaction H36I1 . 
a conventional meson-exchange model, are presented, too. The paper ends with a short summary. In the appendix, we 
provide expressions for the two-boson exchange potential. Furthermore we summarize SU(3) breaking effects which 
arise at NLO from quark mass insertions in the interaction Lagrangian. 
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Figure 1 : Relevant Feynman diagrams up-to-and-including next-to-leading order. Solid and dashed lines denote octet baryons and pseudoscalar 
mesons, respectively. The square symbolizes a contact vertex with two derivatives. From left to right: LO contact term, one-meson exchange, NLO 
contact term, planar box, crossed box, left triangle, right triangle, football diagram. 



2. Chiral potential at next-to-leading order 



The derivation of the chiral baryon-baryon potentials for the strangeness sector at LO using the Weinberg power 
counting has been outlined in Refs. |20, 43. 441 . The NLO contributions for the A^A^ case are described in detail in 
Ref. ll32r . while the extension to the baryon-baryon case has been worked out in Ref. ll45ll . The LO potential consists 
of four-baryon contact terms without derivatives and of one-pseudoscalar-meson exchanges. At NLO contact terms 
with two derivatives arise, together with loop contributions from (irreducible) two-pseudoscalar-meson exchanges. 
The corresponding Feynman diagrams are shown in Fig. [1] 



2.1. Contact terms 



The spin dependence of the potentials due to leading order contact terms is given by 03211 



y 



(0) 



- Cs +CtO-i ■ 0-2, 



(1) 



where the parameters Cs and Ct are low-energy constants (LLCs) that need to be determined in a fit to data. At 
next-to-leading order the spin- and momentum-dependence of the contact terms reads 



y(2) 



= Ciq^ + Czk^ + (dq^ + C^k^) o-j ■ 0-2 + 2^5(0-1 + 0-2) ■ (q x k) 
-i- C6(q ■ cri)(q ■ 0-2) -I- C7(k ■ o-i)(k ■ 0-2) + -Cg((ri - 0-2) ■ (q x k) . 



(2) 



The transferred and average momentum, q and k, are defined in terms of the final and initial center-of-mass momenta 
of the baryons, p' and p, as q = p' - p and k = (p' -H p)/2. The C, (/ = 1, . . . , 8) are additional LECs depending on 
the considered baryon-baryon channel. When performing a partial wave projection, these terms contribute to the two 
5-wave (So, S 1) potentials, the four P-wave (Pi, Po, Pi, P2) potentials, and the Si- Di and Pi- Pi transition 
potentials in the following way 



VCSo) 



4-7t(Cs - 3Ct) + 7r(4Ci +C2- I2C3 - 3C4 - A-Ce - Ci)(,p^ + p'^) , 

vCSi) = 47r(Cs+Cr) + ^(12Ci+3C2 + 12C3 + 3C4+4C6 + C7)(/+p'') 
= 65, + C^s, (/ + //'), 



V('Pi) 



2n 

y 

2n 



(-4Ci +C2 + I2C3 - 3C4 + 4C6 - Ci)pp' = Cif, pp' 



V('Pi) = —{-ACi+C2-ACi + C^ + 2C5-%C6+2Ci)pp'^C,p^pp', 



4 ^n 

y(Vi-'Pi) - — —Cipp'^c.p^^.p^pp', 



(3) 

(4) 
(5) 

(6) 
(7) 



vCPi-'Pi) = -i_^C8p/ = C.p,_3p,pp', (8) 

Vi'Po) = y(-4Ci+C2-4C3 + C4+4C5 + 12C6-3C7)/?/ = C3p„ /?/?', (9) 

Vi'Pl) = y(-4Ci+C2-4C3 + C4-2C5)p/?' = C3p^pp', (10) 

VC^Di-^^i) = ?^(4C6 + C7)//' = C35,_3o,/7'', (11) 

veSi-'Di) = l_^(4C6 + C7);?' = C3s,_3o,;,2^ (12) 

with /? = Ipl and 79' = |p'|. Note that the term proportional to Cg in Eqs. (|7]i and (|7]i represents an antisymmetric 
spin-orbit force and gives rise to spin singlet-triplet transitions (i.e. 'Pi - ^P\). Such transitions cannot occur in the 
A^A^ interaction, unless isospin symmetry breaking is included, and, therefore, this term is absent in the equations 
given in Ref. [32]. However, in general, this antisymmetric spin-orbit term is allowed. Specifically, it does not break 
SU(3) symmetry. 

Assuming only isospin symmetry, the LECs for each spin-isospin channel of the various BB — > BB interaction 
potentials are independent. When imposing SU(3) flavor symmetry one obtains relations between the LECs and, 
thereby, the number of terms that need to be fitted to data gets reduced. The relevant SU(3) structure for the scattering 
of two octet b^ons follows from the tensor product decomposition 8®8 = l©8fl©8j©10*ffilO©27 (for details 
see Refs. HHj). With that one can express all the Cis„,v, C^s,,v, • ■ -, in Eqs. ^-^{v^NN ^ NN, AN -^ AN, 



AN -> I,N, TN -^ I.N) by coeflicients corresponding to the SU(3) irreducible representations: C\ C^% C^\ C'°', 
C'", C^^. The particular combinations of LECs in the various BB -^ BB channels and for the various partial waves 
are summarized in Tab. [1] For example, for the potential in the ' S partial wave of the AJV — > AJV channel we get 

Va^^a^('5o) = ^ [9C?,; + Cf;^ + (9C?]„ + Cf;;(/ + p'')] . (13) 

Note that Tab. [T]gives the weight factors of the various baryon-baryon channels with respect to the irreducible SU(3) 
representations. In addition, it reflects the constraints from the generalized PauU principle. The interaction in partial 
waves like the -^5 1, ^Di, and 'Pi, which are symmetric with regard to their spin-space component, is given by linear 
combinations of coefficients corresponding to antisymmetric SU(3) representations (C^", C''* , C''*), whereas those 
with antisymmetric spin-space part ('5o, ^Po, ^Pu ^Pi) receive only contributions from symmetric representations 
(C^', C^''). The Cg-term induces transitions between singlet and triplet states in the octet-representation 8a and 8.,, 



respectively Il47n . For a detailed derivation of the SU(3) constraints on the LECs see Ref. 1201 or 1481 . 

Due to the imposed SU(3) constraints at LO there are only five independent LECs for the A^A^ and the YN sectors 
together, as outlined in Ref. [20]. Note that without SU(3) symmetry, there would be twice as many. At NLO SU(3) 
symmetry implies that in case of the A^A^ and YN interactions there are eight new LECs entering the S -waves and S -D 
transitions, respectively, and ten coefficients in the P-waves. Note that the sixth leading-order LEC corresponding to 
the singlet representation (C') is present in the strangeness S — -2 channels with isospin 7 = 14911 and there are four 
more LECs that contribute to the S - —2 sector at NLO. 

2.2. Goldstone boson exchange 

The one- and two-pseudoscalar-meson-exchange potentials follow from the SU(3)-invariant meson-baryon inter- 
action Lagrangian 

£mb = tr (b (i/D^ - Mo) b) - f tr (fi/ys {m^, B)) - | tr (fi/yg [u^,,B]) , (14) 

with Df^iB - d/jB + [F^, B], F^ = ^{u^d^u + ud^u^) and u^ - liu^d^u - ud^u'), and where the trace is taken in flavor 
space. The constant Mq denotes the baryon mass in the three-flavor chiral limit. The coupling constants F and D 
satisfy the relation F + D - gA - 1 .26, where gA is the axial-vector strength measured in neutron y6-decay. For the 



Table 1 : SU{3) relations for the various contact potentials in the isospin basis. CV etc. refers to the corresponding irreducible SU(3) representation 
for a particular partial wave f . The actual potential still needs to be multiplied by pertinent powers of the momenta p and p' . 
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ViO 








f = ^So,^Po,^Pu^P2 


^= 'Su'Si-^Du'Pi 


^= 'Pi-^i 
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NN ^NN 
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^10* 


- 




NN^NN 
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cf 
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5 =-1 


AN ^ AN 


1 
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i^(9Cf + C«') 


^{cf.cf) 


-1 r^A 
Vl) ^ 




AN ^I.N 


1 
2 


U-cf-cf) 


\{-cf.cf) 


3 ^8,8,, 




I.N -^ AN 








-I --8,8„ 
V20^f 




I,N ^l.N 


1 
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Ucf^9Cf) 


iicf.cf) 


3 ^8,,8„ 




I.N -^I.N 


3 
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Cf 


^10 


- 



pseudoscalar mesons and octet baryons, collected in traceless 3x3 matrices, 



P = 



(J^ + JL 


n+ 




K^ \ 


V2 V6 








n^ 


-^ + 


JL 


K^ 




V2 


V6 




K 


f 




_2l 








W 



B = 



/-E" A 
V2 "^ V6 


S+ 


p 


2- 


V2 ^ V6 


n 


— u^ 


h" 


2A 



(15) 



V6^ 



we use the usual non-linear realization of chiral symmetry with f/(x) - u^(x) - exp (i y/2P(x) / fo), and /o is the 
Goldstone boson decay constant in the chiral limit. These fields transform under the chiral group SU(3)l x SU(3)r as 
U -^ RUV and B -^ KBK' with L e SU(3)l ,/? e SU(3)r and the SU(3) valued compensator field K = K(L,R, U), 
cf. Ref. 1 50]. After an expansion of the interaction Lagrangian in powers of f one obtains from the terms proportional 
to D and F the pseudovector coupling term 



Xi 



V2 
'2/o 



tr [DByys [d^P, b) + FB^y, [d^,P b]) 



(16) 



which leads to a vertex between two baryons and one meson. In the same way, the term involving the chiral connection 
r^ gives 

^2 = 772tr(iB/[[p,5,P],B]), (17) 

which describes a (Weinberg-Tomozawa) vertex between two baryons and two mesons. 

When writing the pseudovector interaction Lagrangian Xi explicitly in the isospin basis, one gets 

-£i = -fMNnNy^JiTN ■ d^n + iMntyys XH-d^TT 

-U^n [A/rsS + i:/r5A] ■ d^,n - /he^E/tstH ■ d^,n 
-fANK [Ny'ysAd^K + h.c] - f^AK [ErysAd.K + h.c] 
-fmK [N^ysTd^K ■ E + h.c] - f^^K [h/tstS^^ ■ E + h.c] 

-/ee,.S ■ yyslld.Ji - f..,,Z-fys^d,ii . (18) 

Here, we have introduced the isospin doublets 



A^: 



K- 
5 



^0 1 , ^ ^ 



-K 



(19) 



Table 2: Isospin factors I for the various one-pseudoscalar-meson exchanges. 





Channel 


Isospin 


;r 
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?? 


S =0 


NN ^NN 





-3 





1 
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1 
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5 =-1 


AN ^ AN 


1 

2 





1 


1 




AN ^I.N 
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-V3 


-V3 







TN -^I.N 


2 


-2 


-1 


1 






3 

2 


1 


2 


1 



The signs have been chosen according to the conventions of Ref. 114611 . such that the inner product of the isovector S 
(or n) defined in spherical components reads 



S ■ E = ^(-1)"%„S_„, = Z+Z" + E^S" + S"2+ 



(20) 



Since the ori gina l interaction Lagrangian in Eq. (fT&t is SU(3)-invariant, the various coupling constants are related 
to each other by 04611 



fNNn - f, 

/eh.= -(l-2a)f, /he,, 



fNNm^ ^(4q'-1)/, /ank^ -^(l+2a)f. 



/ai.= :^(1-«)/' /22„. = :^(l-«)/' 



-^(l+2a)/, /ea^= ^(4a-l)/, 



f^NK- (l-2a)/. 



(21) 



/S2;r = 2af, 



/aa,,8 = --^(1-a)/, /he^ = -/■ 



Evidently, all coupling constants are given in terms of f = gA/2fo and the ratio a = F/{F + D). 

The expression for the one-pseudoscalar-meson exchange potential is similar to the standard one-pion-exchange 
potential 13211 



V' 



OBE 

B\B2 — >BtiB4 



-JBiB^PjB.BiP Z 5 -Ib,B.^B^B4 ■ 

q^ + nip 



(22) 



Here, nip is the mass of the exchanged pseudoscalar meson. In the present calculation we use the physical masses 
m„, niK, in,j in Eq. (l22l l. Thus, the explicit SU(3) breaking reflected in the mass splitting between the pseudoscalar 
mesons is taken into account. The 77 meson is identified with the octet-state rj^. The isospin factors iBtBi^BjBi are 
given in Tab. |2] 

The two-pseudoscalar-meson exchange potential, built up by a set of one-loop diagrams, is described in detail in 
Appendix A. Relativistic corrections to the one-meson exchange potential that arise at NLO due to differences of the 
baryon masses are discussed in Appendix B. 

2.3. Scattering equation 

In the actual calculation a partial-wave projection of the interaction potentials is performed, as described in detail 
in Ref. OOll . The reaction amplitudes are obtained from the solution of a coupled-channel Lippmann-Schwinger (LS) 
equation: 

T^::r\p",p'^ Vs) = v:::'/(p",p') + ;^ r ^ v^:rv,P) 2 ^^L- C(/^^p; v^) • (23) 

-^ Jo (2ny q^-p^ + 11] "^ 

Here, the label v indicates the particle channels and the label p the partial wave, ^y is the pertinent reduced baryon 
mass. The on-shell momentum qy in the intermediate state, is determined by V^ = ^i^'B "•"??+ v^Bi "^ ^v- 
Relativistic kinematics is used for relating the laboratory momentum pi^b of the hyperons to the center-of-mass mo- 
mentum. 



We solve the LS equation in the particle basis, in order to incorporate the correct physical thresholds. Depending 
on the total charge, up to three baryon-baryon channels can couple. The Coulomb interaction is taken into account 
appropriately via the Vincent-Phatak method 1.51.1 . The potentials in the LS equation are cut off with a regulator 
function, /r(A) = exp -(/5"^ + Pj /^'^U in order to remove high-energy components [28]. We consider cutoff values 
in the range A - 450 - 700 Me V, similar to what was used for chiral A^A^ potentials |28], but anticipate here already 
that the best results are achieved for cutoffs located in the interval 500 - 650 MeV. 

3. Fitting procedure 

For the fitting procedure we consider the same "standard" set of 36 YN data points that have been used in our 



1_P 



previous works Il20t 13611 as also done by the Nijmegen group in their investigations Il37ll . This data set consists of low- 



energy total cross sections for the reactions: Ap -^ Kp from Ref. 15211 (6 data points) and Ref. [53] (6 data points 



irp -^ An [54] (6 data points), I.-p -^ iPn iSM (6 data points), Yp -^ iTp [55] (7 data points), Z+/7 -^ 2+/? [55| 



(4 data points), and the inelastic capture ratio at rest 115 6[ 15711 . Besides these YN data the empirical binding energy of 



the hypertriton ^H is used as a further constraint. Otherwise it would not be possible to fix the relative strength of the 
spin-singlet and spin-triplet S -wave contributions to the Kp interaction. 

We recall that there are in total five independent LECs at LO, that describe the A^A^ and YN interactions, see Tab.[T] 
In Ref. [[200 ^ fit to the FA^ scattering data at LO was presented utilizing these five contact terms. It turned out that 
already in that scenario a fairly reasonable description of the 36 low-energy YN scattering data could be achieved for 
cutoffs A = 550-700 MeV andfornatural values of the LECs. At NLO there are eight new contact terms contributing 
to the S -waves and the ^S \ - ^D\ transition, and ten in the P-waves. As described in Sect. 12. II we impose SU(3) flavor 
symmetry in order to reduce the number of LECs that need to be determined. Without implementing this constraint 
there would be 20 independent contact terms for the AA^ and Y.N systems in the S -waves (and the S -D transitions) 
alone, and, given the low number of data points together with their large error bars, it is simply impossible to fix all 
those LECs by a fit to the available empirical information. 

In the actual fitting procedure, first, we have considered only the 13 LECs in the 5 -waves and the S-D transitions. 
After all, the available YN data consist practically only of total cross sections at low energies and these are predom- 
inantly determined by the 5 -wave amplitudes. The fits are performed for fixed values of the cutoff scale where we 
started with A = 600 MeV. The subsequent fits for other cutoffs were done under the constraint that the results should 
stay as close as possible to those obtained with A = 600 MeV, for the singlet and triplet cross sections separately. This 
procedure is demanded by our requirement to reproduce the hypertriton as mentioned above. 

Note that for the l^p -^ Y^p and E p — > Z p channels the experimental total cross sections were obtained by 



incomplete angular coverage 05511 : 






do-(0) 

o- = ^-!-dcos0. (24) 

— /I — n dcosO 



Following Ref. 113711 . we use cos 6j„i„ = -0.5 and cos fl^ax = 0.5 in our calculations for the H^p — > l,^p and ^.^p -^ I.p 
cross sections, in order to stay as close as possible to the experimental procedure. The total cross sections for the other 
channels are evaluated by simply integrating the differential cross sections over the whole angular region. 
For the capture ratio at rest, r«, we follow the definition of Ref. ll58ll : 

_ 1 crA^-p ^ iPn) ^ 3 (Tr(I.~p ^ iPn) 

''^ 4 o-.CL-p -^ An) + o-,(Y-p -^ IPn) ^ 4 cr,(E-p -^ An) + cr,{Y.-p -^ iPn) ' 

where cts is the total reaction cross section in the singlet '5o partial wave, and cr, the total reaction cross section in 
the triplet-coupled ^S \-^Di partial waves. The cross sections are the ones at zero momentum, but following common 
practice 13711 we evaluate the cross sections at a small non-zero momentum, namely p\.^b - lOMeV/c. 

While the;i'^ fit to the 36 data points allowed us to fix the majority of the S -wave LECs it turned out that concerning 
the ^S 1 partial wave in the 1-3/2 UN channel, solutions with either a positive (attractive) or a negative (repulsive) 
phase shift are possible. This can be understood from the SU(3) structure as given in Table [T]which shows that this 
partial wave is controlled by the "isolated" 10 representation such that the corresponding LECs do not enter in any 



Table 3: The YN contact terms for the 'So and S {- D\ partial waves for various cut-offs. The values of the C's are in 10 GeV the ones of the 
C's in 10'' GeV"'*; the values of A in MeV. 
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450 


500 


550 


600 


650 


700 


'So 


'So 


-0.0893 


-0.0672 


0.00648 


0.1876 


0.6140 


1.145 




'So 


0.2000 


0.1970 


0.1930 


0.1742 


0.1670 


0.1730 




(-27 
'So 


1.500 


1.800 


2.010 


2.200 


2.400 


2.410 




^'S„ 


-0.200 


-0.200 


-0.206 


-0.0816 


-0.0597 


0.1000 


'Sr-'D, 




0.104 


0.541 


0.149 


0.344 


0.499 


0.560 






0.171 


0.209 


0.635 


1.420 


2.200 


2.960 






0.0218 


0.00715 


-0.0143 


-0.0276 


-0.0269 


0.00173 




,-10 


2.240 


2.310 


2.450 


2.740 


2.530 


2.030 






0.310 


0.143 


0.741 


1.090 


1.150 


1.120 






0.373 


0.469 


0.627 


0.775 


0.854 


0.964 




,-10 
^Si-'Di 


-0.360 


-0.429 


-0.428 


-0.191 


-0.191 


-0.122 




/"■lO" 
^Si-'Di 


-0.300 


-0.300 


-0.356 


-0.380 


-0.380 


-0.228 




(-8» 


0.0356 


0.0475 


0.0453 


-0.00621 


-0.00621 


-0.0497 



of the other YN channels. We adopt here the repulsive solution in accordance with evidence from recently measured 
{n^,K^) inclusive spectra related to E -formation in heavy nuclei i59ll60ll6lll . which suggest a repulsive S-nucleus 
single-particle potential |62, 63]. 

We should also mention that we observe some correlations between the values of the S -wave LECs at LO and 
NLO, i.e. C and C in Eqs. O) and ^. This is a consequence of the fact that the fitted l.N cross sections lie all within a 
rather narrow energy interval near threshold so that there i s only a fairly weak sensitivity to the momentum-dependent 
(p^ + p'^) terms. The Kp cross sections alone, which are known over a larger energy range, are not sufficient for 
separating the strength of C and C. 

The limited number (and quality) of differential cross sections and the complete lack of polarization observables 
makes a determination of the contact terms in the f -waves from YN data practically impossible. Therefore, in this 
case and in line with the power counting we assume strict SU(3) symmetry for the contact terms and use the A^A^ 
f -wave phase shifts as a further constraint. In particular, we fix the LECs C^^ and C'" from fitting to empirical 'Pi, 
^Pq, and ^Pi A^A^-phase shifts H in the region of 25 < Tiab < 50 MeV iH]. The ^Pi partial wave is special. Here a 
NLO calculation with the pertinent LEC determined from the low-energy region yields results at higher energies that 
strongly overestimate the empirical phase shifts, see, e.g. Ref. |28]. Such a LEC would likewise lead to a considerable 
overestimation of the Ap cross section around and above the Y.N threshold. In order to avoid this we fix this specific 
LEC from the A^A^-phase shifts in an energy region corresponding to the ILN threshold, namely Tiab ~ 150 MeV. 

Utilizing the A^A^-phase shifts reduces the number of f -wave contact terms that need to be determined in the FA^ 
sector by roughly a factor two. Here the most important constraint is provided by the Kp cross section above the "LN 
threshold, i.e. at around p\.^b ~ 800 MeV/c, which is roughly 10 mb according to experiments ||65[ 16611 . Agreement 



with these data can be only achieved if the contributions from each of the P-waves (Pi, Pq, P\, P2) is kept small, 
which means in turn that the corresponding phase shifts have to be small. The differential cross section for iTp — > A« 
has been measured at two energies near the Y.N threshold and it is sensitive to the P-waves, too. We used this empirical 
information to fix the remaining six P-wave contact terms. But it should be said that this information is not sufficient 
to pin them down reliably. Note that we set the LEC corresponding to the 'Pi-^Pi transition to zero. 

Besides of contact terms with LECs that need to be determined in a fit to data the potential includes also con- 
tributions from one-meson and two-meson exchanges. The latter do not involve any free parameters. The coupling 
constant / = gAl'^fo is fixed by utilizing the values gA - 1.26 and fo ~ fn - 93 MeV. For the FJiF + D)-ratio we 
adopt the SU(6) value a - 0.4. It is close to the empirical value of a « 0.36 - 0.37, as determined recently in analyses 
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Table 4: The YN contact terms for the P-waves for various cut-offs. The values of the LECs are in 10^ GeV "* ; the values of A in MeV. 



A 


450 


500 


550 


600 


650 


700 


Vo 


(-27 


1.47 


1.49 


1.51 


1.55 


1.60 


1.71 






2.50 


2.50 


2.50 


2.50 


2.50 


2.50 


'Pi 


'Pi 


-0.43 


-0.43 


-0.43 


-0.43 


-0.43 


-0.43 




'P, 


0.65 


0.65 


0.65 


0.65 


0.65 


0.65 


'Pi 




-0.096 


-0.063 


-0.041 


-0.025 


-0.012 


0.000 






1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


'Pi 


'Pi 


0.49 


0.49 


0.49 


0.49 


0.49 


0.49 




'P, 


-0.14 


-0.14 


-0.14 


-0.14 


-0.14 


-0.14 




d" 
^'p, 


-0.65 


-0.60 


-0.58 


-0.56 


-0.54 


-0.52 


'Pi-'Pi 


'-'P,-'P, 





















of hyperon semi-leptonic decay data ll67[ 16811 . 

In the spirit of the imposed SU(3) symmetry we keep all contributions from one- and two-meson exchanges, i.e. 
also those from tttj, rjK, KK exchange. The large mass splitting between the Goldstone bosons induces a sizable 
SU(3) breaking in the actual YN potential that is taken into account in our calculation. There is also an explicit SU(3) 
symmetry breaking in the coupling constants as reflected in the empirical values of the decay constants f„, fx and 
/;, |69] which, in principle, should be taken into account in the NLO calculation. However, we have ignored such 
eflrects in the results reported in the present paper. As a matter of fact, we have explored various scenarios in the 
course of our investigation. In particular, we performed fits based on the empirical decay constants. We explored 
also the situation when only tttt exchange diagrams are kept and all two-meson exchanges involving the heavier K, rj 
mesons are omitted. Furthermore, in Ref. [70] we had presented results based on an incomplete NLO calculation, i.e. 
where only the NLO contact terms were taken into account but no two-meson-exchange contributions. In all these 
cases a comparable description of the YN data could be achieved, i.e. with a x^ within 10-15% of the best values 
achieved. Seemingly, all two-meson exchange effects could be absorbed into the LECs and, moreover, one could still 
maintain SU(3) symmetry for these contact terms. A further uncertainty in our calculation is the value of the 77-baryon 
coupling, since we identified the physical 77 with the octet 77s. In our earlier investigation 112 Oil we varied the 77 coupling 
between zero and its octet value and we found very little influence on the description of the data. Thus, we refrain 
from introducing a singlet coupling in the present study. It is possible though that future calculations of hypernuclei 
based on these chiral EFT interactions could indicate a preference for one or the other scenarios and/or yield evidence 
for the need of an explicit SU(3) breaking in the contact terms. 

Finally, let us mention that we did consider also the S -waves for A^A^ scattering. In particular, we fixed the pertinent 
five LECs from a fit to the np ^So and ^S i-^Di phase shifts and mixing parameter in the energy range Ti„b < 50 MeV, 
independently of the YN interaction. Thereby it turned out that the LECs determined from the A^A^ phase shifts are 
incompatible with those needed for the description of the YN data with regard to the SU(3) symmetry. The most 
obvious case is the '^o partial wave, where SU(3) symmetry implies that Vnn = Vsa? (1=3/2), see Tab. [T] so that the 
(hadronic part of the) interaction in the En and S^p channels is unambiguously fixed once the LECs are determined 
from the np phases. However, with LECs determined from the latter channel a near-threshold bound state is obtained 
in the H^p system and, as a consequence, the empirical E+p — > H^p cross section is grossly overestimated. This 
happens despite of the SU(3) breaking in the interaction that arises from the contributions due to one-meson and two- 
meson exchanges and despite of the additional Coulomb repulsion in the "L^p system. Therefore, we must conclude 
that within the scheme followed here a combined quantitative description of the A^A^ and YN sectors based on strictly 
SU(3) symmetric (LO and NLO) contact terms is not possible. Since at NLO an explicit SU(3) breaking in the LO 
contact terms arises anyway, see Appendix B, we considered also a scenario where we took over the NLO contact 
terms from the fit to the A^A^ phase shifts and we fitted the remaining (LO and NLO) S -wave contact terms to the YN 
data. In this case a description of the AA^ and "LN data is possible, but with a noticeably increased ;t'^. Moreover, we 
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Figure 2: "Total" cross section cr (as defined in Eq. (24j) as a function of piab- Tlie expeiimental cross sections are taken from Refs. l52ll (filled 
circles), 15311 (open squares), li65J (open circles), and |66] (filled squares) (Ap —> Ap), from |54!| (Z^p — » An, Z^p — > iPn) and from 15511 
(E^p —> S^p, Z+p — * S^p). The red/dark band shows the chiral EFT results to NLO for variations of the cutoif in the range A = 500,. . .,650 MeV, 
while the green/light band are results to LO for A = 550,. . .,700 MeV. The dashed curve is the result of the Jiilich '04 meson-exchange potential 

a. 



observe the questionable tendency of the Ap amphtude in the S i partial wave to become rather large for momenta 
above the I.N threshold. Thus, we decided to determine all contact terms in the 5 -waves and the S-D transition from 
a fit to the YN sector alone where it turns out that SU(3) symmetry for the LECs can be preserved. 

The values of the contact terms obtained in the fitting procedure for the various cutoffs are listed in Tables [3]and 



HI 



4. Results and discussion 

The results obtained at NLO are presented in Fig. |2] (red/dark bands), together with those at LO (green/light 
bands). The bands represent the variation of the cross sections based on chiral EFT within the cutoff region of 
A = 500 - 650 MeV. Note that in the LO case variations of A = 550 - 700 MeV were considered iQ. For 
comparison also results for the Jiilich '04 [36] meson-exchange model are shown (dashed lines). 

Obviously, and as expected, the energy dependence exhibited by the data can be significantly better reproduced 
within our NLO calculation. This concerns in particular the "E^p channel. But also for Ap the NLO results are now 
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Figure 3: As in Fig.|2] but now tlie experimental cross sections are taken from Refs. pT] (S p — > An, Z p — > iPn), fTs] (X p — > S p), and [74 



well in line with the data even up to the 1,N threshold. Furthermore, one can see that the dependence on the cutoff 
mass is strongly reduced in the NLO case. 

A quantitative comparison with the experiments is provided in Tab. |5] There we list the obtained overall x^ but 
also separate values for each data set that was included in the fitting procedure. Obviously the best results are achieved 
in the range A = 500 - 650 MeV. Here, in addition, the x^ exhibits also a fairly weak cutoff dependence so that one 
can really speak of a plateau region. For larger cutoff values the x^ increases smoothly while it grows dramatically 
when going to lower values. Therefore, in Fig. |2]and in the figures below we show only results based on variations of 
the cutoff within this plateau region. 

A total ;if^ value of around 16 is quite good. Indeed, the best values achieved with phenomenological models, 
say the Nijmegen NSC97 meson-exchange potentials [3/], lie also in that region. We should add that our additional 
requirements that we want to produce a correctly bound hypertriton and that we want a repulsive "LN interaction in 
the isospin 7 = 3/2 channel leads to a shghtly increased ;^f^. Without those constraints we could achieve values which 
are around 5 % smaller In any case, one has to say that one should not overrate the;^'^. Given that there are only 36 
data points the x^ pei" data point amounts to x 0.5 only - which is somewhat low as compared to what one would 
expect from a set of statistically sound data. As a matter of fact, the biggest single contribution to the x^ comes from 
the I.N charge-exchange reaction, see Tab. |5] and specifically from a single data point near threshold that is far off all 
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Table 5: Comparison between the 36 YN data and the theoretical results for the various cutoffs in terms of the achieved x^- The last column in the 
NLO section, denoted 600* , contains result for an interaction where all two-meson-exchange contributions involving the ;;- and/or K meson have 
been omitted, cf. text. 





data 


NLO 


LO 


A (MeV) 


450 


500 


550 


600 


650 


700 


600* 


600 


Kp — > Kp 


Sechi-Zorn [52j 


2.8 


2.1 


2.0 


1.5 


1.6 


1.7 


1.6 


7.5 




Alexander [531 


4.2 


2.6 


1.6 


2.3 


2.4 


2.5 


2.2 


4.9 


iTp — > An 


Engelmann [541 


4.3 


3.7 


3.9 


4.1 


4.4 


4.5 


4.0 


5.5 


Yrp -^ iPn 


Engelmann [541 


5.6 


6.1 


5.8 


5.8 


5.7 


5.7 


5.9 


7.0 


YTp — » l.^p 


Eisele [55] 


2.0 


2.0 


1.8 


1.9 


1.9 


2.0 


1.9 


2.4 


S+p ^ S+p 


Eisele [55J 


0.5 


0.3 


0.4 


0.5 


0.3 


0.6 


0.9 


0.6 


rR 


[56,57] 


0.3 


0.1 


0.2 


0.1 


0.2 


0.3 


0.1 


0.5 


total x^ 




19.7 


16.8 


15.7 


16.2 


16.6 


17.3 


16.5 


28.3 



other data points, see Fig. |2] When one omits this data point the total x^ would be around 1 and the x^ per data point 
would be =s 0.3. Improvements in the order of 5 % on that level are certainlynot significant. In this context let us 
mention thejjf^ of preliminary results of our study that have been presented in ||7U|72{1 is not yet optimal. Specifically, 
there the description of the l.^p cross section was still inferior 

In Table|5]we include also results of an YN interaction where from all two-meson exchange contributions that arise 
to NLO according to SU(3) symmetry only the nn exchange diagrams were kept. All two-meson exchange diagrams 
involving the heavy mesons ;; and/or K were omitted. We performed an exemplary fit within this scenario for the 
cutoff A = 600 MeV and the corresponding ;t'^ values can be found in the column labeled by 600*. It is obvious that a 
comparable fit to the data can be achieved within such a scenario, too. Finally, in the last column of Table|5j A'^ results 
for our LO interaction from Ref. Il20ll (for A - 600 MeV) are reproduced. Evidently, going to NLO allows to reduce 
the x^ by roughly 50 % ! 

A comparison of our results with integrated cross sections at higher energies is presented in Fig. |3] These data 
were not included in the fitting procedure and, therefore, the shown results are genuine predictions of the chiral EFT 
interaction. One can see that the cross sections achieved at NLO are now closer to those obtained from the Jiilich 
meson-exchange potential than the ones at LO, and to some extent they are also more in line with the data. But 
given the large uncertainties in the experiments, even in the fairly recent measurements of the 2 ;? — » iTp [7^ and 
S^p — > Y.^p [74] reactions, precise conclusions are difficult to draw. 



Differential cross sections are shown in Figs. 2] and |5] and compared with available measurements Il54ll55ll73ll74 , 
7511 . Also these data were not included in the fitting procedure (as far as the LECs in the 5 -waves are concerned). 
However, as already mentioned in the preceding section, the differential cross sections for S /? — > An were considered, 
together with the integrated A/7 cross section around piab ~ 750 - 850 MeV/c, in the "by hand" adjustment of the 
LECs in the f -waves. 

As can be seen from Fig. |4] the prediction of the NLO interaction for E^p — > An at piab = 135 MeV/c agrees 
well with the trend of the data [54]. The amplitude is dominated by the ^S\ — > ^Si and ^D\ — > ^S\ transitions 
so that the resulting angular distribution is rather flat. At LO there are significant P-wave transitions (see the phase 
shifts discussed below) that give rise to an enhancement of the cross section at backward angles. Such large P-waves 
arise at LO solely from the one-pion exchange; there are no contact terms yet in those partial waves that would allow 
one to reduce the P-wave amplitudes as it was possible in the present NLO approach. The resulting cross section at 
Piab = 160 MeV/c is very similar. Here the data seem to indicate a clear enhancement in forward direction. It should 
be said, however, that the experimental values shown in Fig. |4]are not the result of a measurement at the specified 
momenta, but rather an average over different momentum intervals. Specifically, for "L^p — > An the data [54J] are 
averages over the intervals 100 < p^,- < 170 MeV/c and 150 < p^.- < 170 MeV/c, respectively. In view of the large 
error bars we refrain here from averaging our theoretical re suits and, following common practice, present predictions 
at the average value of the momenta. The same is also true for the data from Ref. 15511 which represent averages over 
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Figure 4: Diflferential cross section dajd cos 6 as a function of cos d, where d is the cm. scattering angle, at various values of p\j^ . The experimental 
differential cross sections are taken from i54ll (IT p — » An, YT p -^ iPn) and from I55l1 (£"p — » ^^p, ^* p -^ ^*p)- Same description of curves as 
in Fig.|2] 



150 < /?2- < 170 MeV/c (for l.-p -^ iTp) and 160 < pj.* < 180 MeV/c (for l.+p -^ I.+p), respectively. 

Fig. |4] suggests that a somewhat different cocktail of P-wave contributions, like the one predicted by the Jiilich '04 
interaction (cf. the dashed lines), might be more in line with the experimental data. However, we postpone a thorough 
determination of the LECs in the P-waves to future investigations. Here, as a first step, one would try to connect 
our interaction with the effective YN interactions in a nuclear medium as determined from the analysis of y-ray data 
for A hypernuclei |76] via a G-matrix calculation fi77il . The results of such an analysis could provide additional 
and valuable information on the spin-dependence of the YN force and, specifically, on the spin-spin and spin-orbit 
interaction. Indeed, the spin-orbit splitting of the A single-particle levels in nuclei is experimentally well established 
and very small 178117911 and, therefore, should constitute a stringent constraint on the interaction. We expect that then 



also the LEC for the Pi- Pi transition potential can be fixed, which has been set to zero in the present work. 

The data on differential cross sections at higher energies (cf. Fig. |5) are averages over 400 < p^- < 700 MeV/c, for 
S-p -^ Irp and over 300 < pj.. < 600 MeV/c |75] and 350 < /?2;+ < 750 MeV/c |74], respectively, for S+p -^ S+p. 
Also here the predictions of the YN interactions are for the momenta as specified in Fig. |5] 

Scattering lengths for the Ap and l^p interactions in the '5o and ^Si partial waves are summarized in Tab. |6] 
Furthermore we provide results for the hypertriton binding energy. As already said, this binding energy had to be 
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Figure 5 : Diflferential cross section dcrld cos 6 as a function of cos 6, where 6 is the cm. scattering angle, at various values of piat . The experimental 
differential cross sections are taken from [;73i1 (S^p — > S^p), and from L75i1 (filled circles) and 174.1 (open circles) (Z^p — » S+p). Same description 
of curves as in Fig.|2] 



taken as additional constraint in the fitting procedure because otherwise it would have not been possible to fix the 
relative strength of the (5 -wave) singlet- and triplet contributions to the Kp interaction. Tab. |6] lists also results for 
two meson-exchange potentials, namely of the Jiilich '04 model 1.36.1 and the Nijmegen NSC97f potential 1.37.1 . which 
both reproduce the hypertriton binding energy correctly. 

The S^p scattering length in the ^S \ partial wave is positive, as it was already the case for our LO potential, 
indicating a repulsive interaction in this channel. 

Table 6: The YN singlet (s) and triplet (t) scattering lengths (in fm) and the hyp ertriton binding energy, Eg (in MeV). The binding energies for the 
hypertriton (last row) are calculated using the Idaho-N3LO NN potential f27i|. The experimental value for the ^H binding energy is -2.354(50) 
MeV. 









NLO 






LO 


JuHch '04 


NSC97f 
















[20] 


[36] 


[37] 


A [MeV] 


450 


500 


550 


600 


650 


700 


600 




An 


-2.90 


-2.91 


-2.91 


-2.91 


-2.90 


-2.90 


-1.91 


-2.56 


-2.60 


An 


2.64 


2.86 


2.84 


2.78 


2.65 


2.56 


1.40 


2.74 


3.05 


a^" 


-1.70 


-1.61 


-1.52 


-1.54 


-1.51 


-1.48 


-1.23 


-1.67 


-1.72 


.f 


3.44 


3.05 


2.83 


2.72 


2.64 


2.62 


2.13 


2.93 


3.32 


a^'^ 


-3.58 


-3.59 


-3.60 


-3.56 


-3.46 


-3.49 


-2.32 


-3.60 


-4.35 


rf^ 


3.49 


3.59 


3.56 


3.54 


3.53 


3.45 


3.60 


3.24 


3.16 


fl, 


0.48 


0.49 


0.49 


0.49 


0.48 


0.49 


0.65 


0.31 


-0.25 




-4.98 


-5.18 


-5.03 


-5.08 


-5.41 


-5.18 


-2.78 


-12.2 


-28.9 


(IH) Eb 


-2.39 


-2.33 


-2.30 


-2.30 


-2.30 


-2.32 


-2.34 


-2.27 


-2.30 



The Ap scattering lengths predicted at NLO turn out to be significantly larger than those obtained at LO - as 
example for the latter we included the result for the cutoff" A - 600 MeV in Tab.|6] In case of the '^o channel, they 
are even somewhat larger than the values of the meson-exchange potentials. We want to remind the reader that the 
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Figure 6: The Ap 5o and Pi phase shifts 5 as a function of piab- The red/dark band shows the chiral EFT results to NLO for variations of the 
cutoif in the range A = 500, . . . ,650 MeV, while the green/light band are results to LO for A = 550, . . . ,700 MeV. The dashed curve is the result of 
the Jiilich '04 meson-exchange potential l3dl . 



hypertriton binding energy is much more sensitive to the AA^ '^o strength than to that of the ^S i-^Di partial wave, 
as is known from studies in the past llSOilSlll . Thus, the value of the '5o scattering length is strongly influenced by 
our demand to reproduce a correctly bound hypertriton. Interestingly, in the incomplete NLO calculation (i.e. without 
two-meson-exchange contributions) presented in Ref. IITOll a bound hypertriton was achieved with Ap '5o scattering 
lengths around -2.6 fm, i.e. close to the values of the two meson-exchange potentials, cf. Tab.|6] 

The hypertriton results discussed above were all obtained without an explicit three-body force (3BF). It has been 
argued that the variation of the three-baryon binding energy with the cutoff A could provide a measure for the size of 
the 3BF 1 82]. If so one would expect its effect to be somewhere in the range of 10-90 keV, based on the values listed 
in Tab. |6] Formally the first non-vanishing contributions to the 3BF appear at next-to-next-to-leading order (NNLO) 
in the scheme that we follow 12911 . But we want to point out that our present calculation includes already some 3BF 
effects. These are generated automatically in the employed coupled-channel AN-'LN formalism and occur in the 
form of the transition of the A to the E in the intermediate (YNN) state. However, these contributions are two-body 
reducible and, therefore, do not constitute a genuine (irreducible) 3BF. Note that discussions of 3BF eff'ects in the 
strangeness sector in the literature 11231 l24l l25l |261 are often related to the case of an intermediate S, sometimes even 
exclusively. One should distinguish its role from that of an irreducible 3BF which would be generated, for example, 
by the excitation of the S(1385) resonance in the intermediate state - analogous to the 3BF that arises in the standard 
three-nucleon problem due to the A(1232) excitation. 

Calculations for the four-body hypemuclei 1, H and 1, He based on the preliminary version of the NLO interaction 



presented in ||71|] can be found in Ref. Il82ll . That interaction reproduces qualitatively the A separation energies for 



^H and, in particular, it yields the correct ordering of the 0^ and 1^ states. However, a quantitative agreement with 
the experimental information is not achieved. Corresponding computations for the EFT interactions discussed in the 
present paper are in progress Il83ll . 

Finally, let us present predictions for a selection of Ap and I,^p phase shifts, evaluated in the isospin basis. They 
can be found in Figs.|6]-|9] The behavior of the '5o phase shift in the A» channel predicted at NLO is similar to the 
one of the Jiilich '04 model and other meson-exchange potentials ll37ll38ll though may be slightly more repulsive for 
higher momenta, cf. Fig.|6] The 'Pi phase shift is also similar to the result of the Jiilich model and has opposite sign 
as compared to the LO result. Note that this partial wave is the only P-wave where we observed a noticeable cutoff' 
dependence of the results and we counterbalanced this via a smooth variation of the LEC C,° , see Tab.|4] In all other 

Q 

P-waves the value of the additional LEC not determined from the A^A^ sector (C,') is fixed independently of the cutoff 
A. 
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Figure 7: The A/) phase shifts for the coupled S\- D\ partial wave as a function of piab- Same description of curves as in Fig.|6] 



Phase shifts for the coupled ^S \-^D\ partial waves and the mixing parameter e\ can be found in Fig.|2l Evidently, 
the ^S 1 phase shift based on the NLO interaction passes through 90° slightly below the 'LN threshold. However, also 
in the ^D\ phase shift and the mixing parameter we observe an appreciable increase near that threshold. A rise of the 
^5 1 (or ^Di) phase shift beyond 90° is typical for the presence of an unstable bound state in the "LN system 11841851] . 
see also the discussion in [86]. In case of the l^A^ interaction at LO and the Jiilich '04 model [36] none of the phases 
pass through 90° and an ordinary cusp is predicted. Such a behavior is caused by an inelastic virtual state in the 
'L'N system. It should be said, however, that the majority of the meson-exchange potentials OSl l37l 1381] produce an 
unstable bound state, similar to our NLO interaction. The only characteristic difference of the chiral EFT interactions 
to the meson-exchange potentials might be the mixing parameter e\ which is fairly large in the former case and close 
to 45° at the "LN threshold, see Fig. Q It is a reflection of the fact that the pertinent Ap T-matrices (for the ^S\^S\, 
^D\^D\, and ^S\'^^D\ transitions) are all of the same magnitude. Indeed, these amplitudes yield very similar 
contributions to the Kp cross section in the vicinity of the SA^ threshold. 

In this context let us mention a recent experimental paper where the energy region around the Y.N threshold was 
investigated in the reaction pp — > K^hp via a measurement of the Kp invariant mass 11871] and where a pronounced 
structure was observed. For a discussion and summary of older measurements providing evidence for a strong en- 
hancement of the Ap amplitude near the "LN threshold see Ref. (Saj. 

Predictions for the ^P partial waves of the Kp system are displayed in Fig. |8] One can see that the ^Pq and ^P\ 
phase shifts are reduced at NLO as compared to those obtained atLO while they are larger in case of the ^P2- Note that 
the behavior of the NLO results is strongly influenced by the LECs as fixed from the corresponding A^A^ partial waves 
because, according to SU(3) symmetry, the pertinent coefficient (C^^) dominates also the AA^ — > AA^ interaction, 
see Tab. [T] Obviously, there are sizable quantitative differences between the results for the EFT interaction and the 
meson-exchange potential. 

Results for the S^p system are shown in Fig. |9] where we restrict ourselves to the 5 -waves. We have switched 
off the Coulomb interaction for the computation of the phase shifts so that the presented results are those for the YN 
I - ijl channel. There is no coupling to the AN system and therefore the phase shifts are elastic in the momentum 
region considered. 

Both partial waves are quite interesting. First, with regard to the '5o, strict SU(3) symmetry implies that Vnn = 
V-E,N, see Tab. [T] so that in an exactly SU(3) symmetric world the corresponding pp and Y.^p phase shifts would be 
the same. In our calculation we break the symmetry already on the potential level by using the physical masses of the 
exchanged pseudoscalar mesons and in addition by using the physical baryon masses when solving the LS equation 
(l23l l. As already mentioned in Sect. [3] we had to introduce also an SU(3) breaking into the contact terms. It turned 
out to be impossible to describe the pp '5o phase shifts and the S^p cross sections with a consistent set of LECs that 
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Figure 8: The Kp P-wave phase shifts 5 as a function of piab- Same desciiption of cui'ves as in Fig.|6] 



fulfill SU(3) symmetry, at least on the level of our NLO calculation. As a matter of fact, the '5o amplitude of our 
NLO interaction alone saturates already more or less the experimental E^p cross sections. Employing the LECs as 
fixed from a fit to the pp phase shifts yields a potential that is much more attractive and that produces a near-threshold 
bound state in the S^p system and, consequently, cross sections that are roughly four times larger than experiment. 

Apparently meson-exchange interactions like the Nijmegen potentials are able to describe the A^A^ and YN systems 
simultaneously, without any major obvious SU(3) breaking. However, usually a special (phenomenological) treatment 
of the short-r anged p art of the potential is required, as discussed, for example, in Ref. OSll . In the potentials of the 
Jiilich group 135113611 SU(3) symmetry is broken via the employed vertex form factors. 



The coupled S \- Di partial wave of the l^p system has a strong influence on the properties of the E in nuclear 



matter Specifically, a repulsive S single-particle potential in nuclear matter 11621 16311 . as supported by present days 
experimental evidence JSOlroOlrolll . can only be achieved with a repulsive ^S\ interaction in the I - 312 channel. 
In the course of our investigation we found that we can fit the available YN data equally well with an attractive or a 
repulsive ^S \ interaction. The difference in the achieved ;\f^ is marginal as already pointed out above. In view of the 
SU(3) structure given in Tab. [T]this may be not too surprising. The ^^i partial wave of the 1,N I = 3/2 channel resides 
in the 10 representation which does not contribute to any of the other A^A^ and YN systems. Of course, its contribution 
enters indirectly because the measured (physical) reactions YTp — > iTp and YTp — > iPn involve amplitudes that result 
from combinations of the Y,N I - ill and 1-1/2 interaction potentials. 

Our NLO interaction produces a moderately repulsive ^S i phase shift as can be seen in Fig. |9] comparable to the 
one predicted by the LO potential. For the latter, calculations of the 2 single-particle potential have been performed 
1 8811 and indicate a values of U^(k = 0) ~ 12 MeV at nuclear matter saturation density. 

Recent lattice QCD calculations 1 19] suggest a much more strongly repulsive ^S i phase shift in the 1,N 1-3/2 
channel, when extrapolating the lattice results obtained for m„ ^ 389 MeV to the physical pion mass. But within 
our framework we cannot accommodate a much more repulsive ^S i amplitude. Any sizable increase in the repulsion 
would yield a ^S \ amplitude which practically saturates the experimental l,^p cross section alone and, consequently, 
there would be no room anymore for the contribution from the spin-singlet amplitude - which is likewise large as 
discussed above. Thus, a more strongly repulsive ^S i phase shift would immediately result in a dramatic deterioration 
of the achieved ;^f^. 



5. Summary and outlook 



Chiral efifective field theory, successfully applied in Refs. 11271 12811 to the A^A^ interaction, also works well for the 
baryon-baryon interactions in the strangeness S - -1 (AN - 'LN) and S - -2 (AA - SA^ - YZ) 1201 14911 sectors. 
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Figure 9: The Z^p S-wave phase shifts 5 as a function of piab- Same description of curves as in Fig.|6] 



As shown in our earlier work Il20ll . already at leading order the bulk properties of the AN and 'LN systems can be 
reasonably well accounted for The new results for the YN interaction presented here, obtained to next-to-leading order 
in the Weinberg counting, are very encouraging. First there is a visible improvement in the quantitative reproduction 
of the available data on AjV and I,N scattering and, secondly, the dependence on the regularization scheme is strongly 
reduced as compared to the leading-order result. Indeed the description of the YN system achieved at NLO is now on 
the same level of quality as the one by the most advanced meson-exchange YN interactions. 

At the considered order there are contributions from one- and two-pseudoscalar-meson exchange diagrams and 
from four-baryon contact terms without and with two derivatives. SU(3) flavor symmetry is used as guiding principle 
in the derivation of the interaction. This means that all the coupling constants at the various baryon-baryon-meson 
vertices are fixed from SU(3) symmetry and the symmetry is also exploited to derive relations between the contact 
terms. Furthermore, contributions from all mesons of the pseudoscalar octet (n, K, iff are taken into account. The 
SU(3) symmetry is, however, broken by the masses of the pseudoscalar mesons and the baryons for which we take 
the known physical values. 

Given the presently available data base with its still large uncertainties, we are able to achieve a combined de- 
scription of the AA^ and ILN systems without any explicit SU(3) breaking in the contact interactions. However, we 
found that a simultaneous description of the A^A^ interaction with contact terms fulfilling strict SU(3) symmetry is not 
possible. Here the strength of the contact interaction in the 27-representation that is needed to reproduce the pp (or 
np) '^o phase shifts is simply not compatible with what is required for the description of the empirical S^/? cross 
section. 

In any case, it is likely that future (and more precise) data will demand to depart from SU(3) symmetry in the 
contact terms even with regard to the AN and ILN interactions. Especially studies of few- and many-body systems 
involving hyperons, which can be done in a consistent way in the framework followed in the present work, could 
provide evidence for the need of an explicit SU(3) breaking. So far reliable microscopic calculations that utilize 
directly the elementary YN interaction are only possible (and have been done) for systems with at most four baryons, 
namely within the Faddeev- Yakubovsky approach 18911 . However, it is expected that new approaches that have been 
developed and refined over the past few years and that are successfully applied in studies of ordinary nuclei allow one 
to study also hypemuclei with a larger number of baryons with comparable accuracy. Thus, we consider the present 
investigation as a first and exploratory step towards a more thorough understanding of the baryon-baryon interaction 
in the strangeness sector 
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Appendix A. Two-pseudoscalar-meson exchange contributions 

In this section we present the next-to-leading order contributions from two-pseudoscalar-meson exchange shown 
in Fig. [1] The calculation of these potentials was done according to the rules of SU(3) heavy-baryon chiral perturba- 
tion theory in the center-of-mass frame and in the isospin limit. Ultraviolet divergences are treated by dimensional 
regularization, which introduces a scale A. These divergences are parametrized in an R-teixn which is absorbed by 
contact terms. In the used renormalization scheme it is defined as 



R = 



d-A 



+ Te - 1 - In (4;r) , 



(A.l) 



with the space-time dimension d. 

As for the one-pseudoscalar-meson exchange, Eq. (l22t . the two-pseudoscalar-meson exchange potentials are given 
by a general expression, where the proper meson masses have to be inserted, and which has to be multiplied with 
appropriate SU(3) factors A^. We display this factor next to the Feynman diagram and in the corresponding tables. The 
factors contain coupling constants and isospin factors and are different for each combination of baryons and mesons. 

In the following we will show the results for the five diagram types one after another 

Appendix A.l. Planar box 
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Figure A. 1 : Planar box 

The planar box. Fig. lA.ll has an irreducible part and a reducible part coming from iterating the one-meson exchange 
to second order The reducible part is generated by solving the Lippmann-Schwinger equation and, therefore, is not 
part of the potential. The irreducible part of this diagram can be obtained by regarding the residues of poles of the 
meson propagators, but leaving out poles of the baryon propagators. One obtains a central potential ( IVc), a spin-spin 
potential {(T\ ■ (T2V5) and a tensor-type potential {<T\ ■ qcr2 ■ qVr). With the momentum transfer q - |p' - p | and the 
masses of the two exchanged mesons, mi and m2, the irreducible potentials can be written in closed analytical form. 



T zPlanar box , ^ 



A^ 



1927r2 



[m\ - ml) 



+ 16 {t7i\ + ml) + {23q^ + 45 [ml + ml)) iR + 2ln 



■\Jm1m2 



2 2 



(uqU U - mlf - 9q^ U + m^)) In ^ + -^i22>q' - fcL^ 



+ 56 [ml + ml) q^ + 8 \ ^ [ml -ml) +2 [iW^ + 22m?m^ + 21m^) L (q) 



(A.2) 
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Table A. 1: Isospin factors I for planar box diagrams. BuBjr indicates the two baryons in the intermediate state and ttjt etc. the exchanged pair of 
mesons M1M2, cf Fig. lA.ll 
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(A.3) 



(A.4) 



The relation {(T[ x q) ■ (cr2 x q) = q ctj • cr2 - (cri ■ q) (cr2 ■ q) is exploited for the connection between the spin-spin 
and tensor-type potential. The isospin factors I can be found in Tab. lA.ll Two-meson exchange contributions that 
involve a single K (or K) lead to an interchange of the nucleon and the hyperon in the final state. The recoupling of 
the corresponding isospin states yields a factor (-1) for some of the transitions that is already included in the values 
given in Tab. lA.ll The same applies to the Tables given below. In this context let us mention that for diagrams with 
an interchange of the nucleon and the hyperon in the final state, likewise an appropriate treatment of the spin-space 
part is required. In particular, the momentum transfer is then given by ^ = |p' -H p |. 

Note that the potential given above and also the following potentials are finite for ^ — > 0. Terms proportional 
to 1/^2 and/or l/q'^ in Eqs. ( IA.2b and ( IA.3b are cancelled by corresponding terms in the functions L{q) and w(q) of 
Eq. (IA.4b in the hmit of small q. We perform an expansion of the potentials in a power series for small q so that these 
cancellations can be taken into account analytically and we obtain stable results in the numerical calculations. Fo r 
equal meson masses some terms in the potentials vanish and the expressions reduce to the results in Refs. \\32 
This is the case in the A^A^ interactions of Refs. 



m 



32L 1271 12811 based on chiral EFT, where only contributions from 



two-pion exchange are taken into account. 

In the actual calculations only the non-polynomial part of Eqs. ( IA.2I ) and ( IA.3b is taken into account, i.e. the pieces 
proportional to L{q) and to 1 /q^ and 1 /q*. The polynomial part only renormalizes the LO and NLO contact terms and, 
therefore, is not considered. The contributions involving the regularization scheme (i.e. that depend on R) are likewise 
omitted. As already said, their effect is assumed to be also absorbed by the contact terms and a renormalization of the 
coupling constants, see, e.g., the corresponding discussion in Appendix A of |32] for the A^A^ case. We want to remark 
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that the majority of those terms omitted involve the masses of the pseudoscalar mesons and, therefore, generate an 
SU(3) symmetry breaking. Thus, the SU(3) symmetry imposed on our contact interaction (at least for AA^ and I,N) is 
understood as one that is fulfilled on the level of the renormalized coupling constants. 

All statements above apply also to the other contributions to the potential described below. 
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Figure A. 2: Crossed box 
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The crossed box diagrams. Fig. IA.2J yield a central, a spin-spin, and a tensor-type potential. Due to the similar 
structure but different kinematics, the potentials resulting from the crossed boxes are the same as those of the planar 
box, up to a sign: 



ycros,sedbox(^) ^ _ ^planar box ^^^ ^ 

ycrossed box^^-, ^ _ j_ ^crossed box^^-, ^ yplanar box^,^^ ^ 
1 



(A.5) 
(A.6) 



Note that there is no iterated part in case of the crossed boxes. The corresponding isospin factors I can be found in 
Tab.|A21 
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Figure A. 4: Right triangle 

N - fB,BiM,fBiB,M2J^B,B2^B,Bi 

The two triangle diagrams. Figs. IA.3l and lA.4l lead to equal potentials, but with different SU(3) factors. They con- 
tribute only to the central potential and the corresponding expression is given by 
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Table A.2; Isospin factors I for crossed box diagrams. B,7B,> indicates the two baryons in the intermediate state and Tin etc. the exchanged pair of 
mesons M1M2, cf. Fig. lA.2l 
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Table A. 3: Isospin factors I for triangle diagrams with the BBMM vertex at the left baryon. B, denotes the baryon in the intermediate state and nn 
etc. the exchanged pair of mesons Mi M2, cf. Fig. lA.3l 
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The isospin factors I can be found in Table lA.5 
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Table A. 4: Isospin factors I for triangle diagrams with the BBMM vertex at the right baryon. B; denotes the baiyon in the intermediate state and 
nn etc. the exchanged pair of mesons M^Mi, cf. Fig. lA.4l 
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Table A. 5: Isospin factors I for football diagrams, nn etc. indicates the exchanged pair of mesons M1M2, cf. Fig. lA.5 
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Appendix B. SU(3) breaking 

Appendix B.l. SU(3) breaking in the contact terms 

In addition to the SU(3) symmetric contact terms given in Sect.|2]that arise at NLO, there are further contact terms 
at this order that lead to an explicit SU(3) symmetry breaking. These terms contain new, i.e. additional, low-energy 
constants. As already mentioned in Sect. |3] the lack of experimental data makes it practically impossible to fix those 
contact terms and, therefore, we decided to set all the corresponding constants to zero. However, for completeness 
and for future reference, we summarize here the structure of the pertinent contributions. 

First there would be, in principle, relativistic corrections (l/Mg) to the leading order contact terms 1I20I1 . 

£1 = C[ tr [iB"{VE),¥(VB)p) , £2 = C^ tr {b"B''(VB)is(VB)„) , £3 = 4 tr (B"(r'B)„) tr {b''(VB)is) , (B.l) 

which break SU(3) symmetry because of different baryon masses. Here a sum over the different elements of the Dirac 
algebra, P e {l,')';i,')'5,y5y;j,cr^v), is implied. The indices a and/? are Dirac indices. However, since the corrections 
to the baryon mass in the chiral limit are of order 0{q^), explicit symmetry breaking due to these corrections does not 
appear up to NLO. 

However, NLO contact terms with an insertion of the external field ;^f, which is of order 0{q^), are possible. In the 
case of baryon-baryon scattering that field amounts to 



X = 27?o 



nu 





0] 




/ 2 








md 





a! 





m; 








m, 











2ml 



(B.2) 



The following baryon-baryon contact terms with insertions of ;^' are possible 114 811 : 

£i^C[tr{B''x(r'B)J''(rB)^), 
£2 = C'2tv{B''(rB)„xB''(rB)f,) , 
£3 = C^tr(B";rfi^(r'B);;(rB)„) + tr(B"B^(r%;ir(r'7J2)„) 
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X4 = C\\x[B"BPx(T'B)p(rB\) , 

Xs = C>(B"F(rB)^(rB)„;(r) , 

Xe = C>(B"(rB)„;^r)tr(B^(rB)^) , 

X7 = C; tr (B";tr ) tr {(rB),B''(rB)^) + tr (B"(rB)„B^) tr {(rB)/,x ) , (B.3) 

and lead to an explicit SU(3) symmetry breaking linear in the quark masses, since we use niu - nid + nis- Using these 
Lagrangians one obtains terms that contribute to the '5o and ^Si partial waves (only). In the following we list the 
results for different transitions (with isospin /). 
NN -^ NN, 7 = 0: 

V('5o) = 0, (B.4) 

VeS 1) = 4n[4Cl^ml + C?^_ (8m| - 4ml) + C%^ (Sm^ - 4ml) \ 

= m^C"'' , (B.5) 

NN ^ NN, 7=1: 

VeSo) = 4n[4C}^ml + C?,^ (8m| - 4ml) + Cf,„ (Sm^ - 4ml) ] 

= mlC^' , (B.6) 

vCSi)^0, (B.7) 

AN -^ AN, I = 1/2: 

Vi'So) =4;r[iC,',^ (4m^ + ml) + C?,_^ (3m| - |mj) + C?^^ (|m^ - ml) + \C%ml 

+ K?,„ (2< - m^) + fCf,^ (5m| - 2^2) + |C/^^ [m\ - ml) ] 

^^ml{9C^^ + C'), (B.8) 

VeS,) =47r[iC3',_ (4m^ - ml) + iC?,_ (7m| - 4ml) + C^s, {"^^ ' ^l) + IC^s^ml 

+ C?^_ (3,4 - im2) + §C3^^_ {5ml - 2ml) + 4C3\_ (m^ - m^) ] 

= im2(c"'' H-C*^") , (B.9) 

AN ^I.N, 1= 1/2: 

y('5o) = 4n[c}^ml + C?^ m^ + C?^^ [ml - 2ml) " K'.o'"^ + K?5„ K " ^'"i) " ^C's. H + '"i) 1 

= ^m2 (C^' - C27) + (^1 _ „,2) Cj , (B.IO) 

y(35 1) = 4n[ - Cl^ml - C^ml + C%^ (2m^ + m^) + fC^, m^ + C35^^ (3m^ - fm^) + 2Cj,^ [ml - ml) ] 

= ^ml[C^°' -C^'') + [ml-ml)C2, (B.ll) 

I.N ^l.N, 7= 1/2: 

y('5o) = 4n[ - Cl^ml - C\ml + 3Cf^ m^ + |Cf^ mj + \C%^ (2m| - m^) + 2C%m'^ 

= jt^m^ (C2^ + 9C**') + (m^ - ml) C, , (B.12) 

y(35i) = 4n[cl^ml + C^^ m| + 3C3\ m^ + ^C^ml + ^C!^^ (2m^ - ml) + 2C%ml\ 

= im2(C'"*+C''") + K-'w')C4, (B.13) 
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I.N -> lA^, / = 3/2: 

VeSo) = ^n[2C}^ml + lC\m\ + 2Cf^ml] 

= m2c27 + (,4-,„2)c5, (B.14) 

Vi^Si) = 4n[ - 2Cl^ml - 2C\m\ + 2&,^m\\ 

^mlC^^. (B.15) 

We introduced here appropriately redefined constants C so that the SU(3) breaking is clearly visible. In case of 
flavor symmetry where m„ - mk the leading order SU(3) relations, cf. Tab. [T] are obtained and the constants can be 
absorbed in the leading order contact terms. If m„ + niK one obtains additional (though suppressed) constants that are 
proportional to m^ - m^. 

Appendix B.2. SU(3) breaking in the OBE contribution 

At NLO and NNLO there are corrections to the one-meson exchange potential due to dififerences in the baryon 
masses. Energy conservation leads to 



^p2 + Ml + ^p2 + Ml = ^p'2 + Ml + ^p'2 + Ml , (B.16) 

and, therefore, in some cases p^ i^ p'^ and/or ^o t^ 0, where 

^" = A£ = El, - El = El - El, . (B. 17) 

Using Mb, - Mq + O \jr\ JSOl] and performing an expansion in 1 /Mo one obtains 

1 



^B,B2^B3B4 ~ ~ fBiBiPfB.BiPJ^BiBi^BiBt 



/7 ? 

P P / _/ _' _ _x , ?0 



o"! ■qo-2q 

(B.18) 



if ir / / / \ ^y> / f / \ 

+ A.. 7 (o-i ■ P o"2 ■ P -o-ip(r2P) + -r-{,(rip(r2P - o"! ■ p 0-2 ■ p ) 

4-Mf, Mo 



'0 



The first term gives rise to the leading order tensor potential, see Eq. (l22l) . but with a shift q^ —> 1^ ~ ?o caused by 
the mass diff'erences of the baryons, i.e. ^" « AM where AM - (Mb, + Ms_, - Mb, - Mb2)/2 ^M- The last two terms 
in Eq. ( IB. 18b give a formal contribution beyond LO. The term proportional to (p'^ ~ P) contributes, in general, only 
off'-shell. An exception are transitions where the baryon masses in the initial state are not equal to those of the final 
state, cf. Eq. (IB. 161 ). For the YN interaction considered here this is only the case for the Van^-ln transition potential. 
In the present study we have neglected all these corrections. 

There are also deviations of the meson-baryon coupling constants from the SU(3) values which, in principle, 
should be taken into account in a NLO calculation. Specifically, there is an explicit SU(3) symmetry breaking in the 
empirical values of the decay constants 16911 . 

A = 92.4MeV, /, = (1.19 ± 0.01)/,, /^ = (1.30 ± 0.05)/, . (B.19) 



A somewhat smaller SU(3) breaking occurs also in the axial coupling constants, see Il67ll68ll91[l but also Il92ll93[l . All 
these eff'ects are likewise not taken into account in the present study. Rather we use the standard SU(3) relations for 
the baryon-baryon-meson coupling constants Eq. (l2ll with the values gA = 1.26 and fa ~ fn - 93 MeV. 
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